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Abstract. The effect of the correlation induced by global momentum conservation on the 
two-particle distribution in nucleus-nucleus collisions is discussed, with a focus on the generic 
case of collisions with a non-vanishing impact parameter. 



1. Introduction 

As is well known, the distribution of particles in the final state of a collision between elementary 
particles results from the interplay of two ingredients |T]. First, there is a "phase space" part, 
which reflects the kinematic constraints arising from energy and momentum conservation. Then, 
there comes the dynamical component, namely the detailed interaction(s) between the "initial- 
state" particles. This is arguably the most interesting part, inasmuch as one does not expect any 
violation of energy-momentum conservation, where the modeling enters — either first-principle 
based modeling for the interaction of two elementary particles, or more effective descriptions as 
required when more particles are involved — together with its set of assumptions and parameters. 
In particular, the aim of experimental investigations is to determine, or at least to constrain, 
those possible models, taking the conservation of energy and momentum as granted. For that 
purpose, it is necessary to have a good control of the "physically trivial" kinematic component. 

For single-particle distributions, the effect of energy- momentum conservation is typically to 
limit the value of the particle momentum and to result in a depletion in the yield close to the 
boundary: the effect amounts to a multiplicative factor. In the case of the joint distributions 
of two, three or more particles, the influence of the phase-space constraint is usually less 
straightforward. If there are only a few final-state particles, the joint distributions are quite 
constrained and the corresponding possible final states are conveniently represented on Dalitz 
plots [Ij. However, in high-energy collisions, in particular of heavy nuclei, where many particles 
are emitted, such an approach becomes impractical. One has instead to estimate the effect of 
total momentum conservation on a statistical basis, be it for the induced correlation between 
two [2j or more than two ^ particles. 

Here, I shall briefly recall how the multiparticle correlation induced by the conservation of 
total momentum can be derived in the limit of a large number of emitted particles [3]. Then 
I shall focus on a case of interest for heavy-ion collision studies, namely when the emission of 
particles in the plane transverse to the beam is anisotropic (i.e., in the presence of so-called 
"anisotropic flow"). 



2. Cumulants from momentum conservation 

Let us consider a collision in which particles are emitted, viewed in the center-of-mass frame 
of the particles, so that the sum of their momenta vanishes + • • • + pat = 0. This relationship 
induces a correlation between the momenta of M particles chosen arbitrarily among the 
ones: technically, the corresponding joint M-particle probability distribution does not factorize 
into the product of the M single-particle probability distributions, but involves non-vanishing 
cumulants [1] of all orders 

The strength of the M-particle correlation depends on both M and the total number of 
particles A^. In the simplest case oi N = 2 final-state particles, the correlation is maximal, since 
P2 = — Pi; the corresponding two-particle probability distribution /(pi,P2) = /(pi)'^(pi +P2) 
is obviously not factorizable. In the large-A^ limit, one can prove with the help of a saddle- 
point approximation that the M-particle cumulant scales like 1/N^^~^ [3]. A convenient and 
systematic approach to show that, as well as to obtain the explicit expressions of the cumulants, 
consists in using a generating function of the joint multiparticle distributions, the logarithm of 
which generates the cumulants. A few steps allow one to write this generating function as the 
integral of the exponential of A^ times some function of the integration variable k. The 

procedure to derive the successive cumulants then consists in computing to a given order in 
powers of 1/A^ the position of the maximum of i.e. the saddle point, then to calculate the 
value of J-" at this maximum so as to perform the saddle-point integration 

Admittedly, the cumulants induced by total momentum conservation are in general small 
when A^ is large. However it is worth keeping their existence in mind, since they can become 
significant in some regions of phase space. At the same time, the correlations induced by 
other more dynamical phenomena which one attempts to investigate might conversely be small, 
and thus not necessarily significantly larger than the "trivial" kinematic ones. Thus, it has been 
argued that some influence of momentum conservation on the measured two-particle short-range 
correlations of identical pions in pp collisions at RHIC energies can be evidenced |6j. In heavy 
ion collisions, where A^ is larger, momentum-conservation induced correlations are even smaller, 
yet they could play a non-negligible role in some studies of small signals. Since it is important 
to have a good idea of what their effect might look like, so as to try to identify and subtract 
similar patterns in the measured correlations, I shall now discuss further these cumulants. 

3. Azimuthally-dependent cumulants and distributions 

The most general expressions of the two- and three-particle cumulants due to the conservation of 
global momentum to leading order in powers of 1/A^, derived according to the method sketched 
above [3], can be found in reference [7] (equations (3.4) and (3.5)). Neglecting the components 
of the momenta along the beam directior|^, I shall focus on the constraint from transverse 
momentum conservation. 

In the heavy- ion context, the mean square momenta along the nucleus-nucleus impact 
parameter (i.e., in the reaction plane) (p^) and perpendicular to it in the transverse plane, 
(Py), are generally unequal: this is the celebrated anisotropic expansion. To account for the 
phenomenon, let me introduce the coefficient V2 = {p^ —Py)/{Px Note that while a 

^ Using probability distributions, instead of the distributions themselves, is much simpler, for it circumvents the 
issue of choosing a "proper" normalization — e.g. of particle-pair yields — when the multiplicity fluctuates 
from event to event (see the discussion in reference [5], section III). Yet, one should not forget that the problem 
is always present in a real experimental correlation measurement. 

^ In the absence of an experimental estimate of the mean square momentum along the beam as compared to along 
the transverse directions, this approximation can tentatively be justified by the reported quasi-independence of 
the average transverse momentum (pr) on rapidity [8]. 

This definition differs from that of the usual elliptic-flow coefficient V2 = {{Px ~ Py)/{Px +Py)), yielding values 
larger by about a factor 2. 



fourth-harmonic modulation V4 of the single-particle distribution has also been evidenced at 
RHIC, I need not consider it here, as it would affect the cumulants induced by momentum 
conservation only from the four-particle cumulant onwards. Introducing V2 into the two-particle 
cumulant expression, one can recast it as 



/c(Pri,PT2) = - j,r/ 2\n 7^ [cOs((/?2 - ^^l) " ^2 COS (991 + ip2 - 2^r)] , (1) 



where $ij is the reaction-plane azimuth. This two-particle cumulant depends not only on the 
relative angle A(pi2 = ip2 — (fii between the particles, but also on the absolute orientation of 
the particle pair with respect to the reaction plane. In turn, the joint probability distribution 
/(Pti,PT2) = f{v'^'^",PTi,PT2,^^i2) depends on the "pair angle" (/jP^"' = {ipi + 992)/2 as well. 
Given equation ([T]), computing the two-particle distribution is straightforward. Keeping only 
the second-harmonic modulation ("elliptic flow") of the single-particle distributions, one finds 

/('^P"^P^l,PT2, AVP12) = ^ [1 + 2^Pf cos2(<^P-'-- ^r) + 2z;Pf sin2((^P-'-- + (2) 

where the Fourier "pair-flow" coefficients [5\ that characterise the azimuthal dependence of 
/(Ptii PT2) depend on pTi, PT2 ^^^^ ^^12- Considering only terms up to 0{v2/N) and ^(fDllI 

<f (PTi,PT2, A9912) ~ [1 - 2t;2(l)t;2(2)cos(4A(/Pi2)] ^(l) + ^^2(2)] cos(2A99i2) + ^§^^2, (3) 

v^f{pTi,PT2,^Vi2) ^ [1 - 2t;2(l)?;2(2)cos(4A(/7i2)] ^(l) - ^^2(2)] sin(2A(/Pi2), (4) 

where ^2(1), "^2(2) are shorthand notations for V2{pti) and V2{pt2)i respectively. There are also 
higher harmonic terms (ff c"^, V4,^"- ■ ■), which are however smaller by at least a factor V2- The 
effect of momentum conservation is actually subleading in the sine term (jl]), which reflects the 
non-invariance of the system under the (p^^^^— —[ip^^^^— symmetry when the particles 

in the pair are different and/or have different transverse momenta [5], and is mostly due to 
anisotropic flow. On the other hand, momentum conservation affects the cosine term ([3]), which 
is non-zero even for particles with vanishing elliptic flow. Since ^2 > at relativistic energies, 
V2^" is positive for small-angle (|A(/3i2| < f ) or large-angle (|A(^i2| > ^) pairs. This means 
that the yield of such pairs is larger in the reaction plane ((p^^^" f» than perpendicular 
to it. Since V2{pt) grows with increasing transverse momentum, as does the second term in 
the right-hand of equation ^2'^''^ increases with both p^i and pT2y the anisotropy in 
the pair yield increases with the particle transverse momenta. Stated differently, for a pair 
of particles close in azimuth {Aipi2 close to 0), the anticorrelation ([1]) is smaller (resp. larger) 
when the pair azimuth is along (resp. perpendicular to) ^r, so that the pair yield is less (resp. 
more) "suppressed" by momentum conservation: there are overall more particles to balance the 
pair momentum along $r than out-of-plane. Conversely, for a pair of back-to-back particles 
(A(/?i2 ~ vr) momentum correlation induces a positive correlation ([1]), which is larger if both 
particles lie along the reaction plane (resulting in ipP^^^ fa ^r ± ^) than perpendicular to it. 

Rephrasing the above in yet another manner, one can investigate the conditional probability 
to find an "associated" particle (transverse momentum PT2) given a "trigger" particle (pTi), by 
dividing the pair distribution /(pti,PT2) by the single-particle distribution f{pTi)- One then 
finds that for a trigger along the reaction plane {tpi ~ $_r), there is a higher probability for 
associated particles close or away (A(/3i2 ~ or vr) in azimuth than around A(pi2 ~ it^. On 
the opposite, if the trigger escapes the system perpendicular to ^r, the conditional probability 



* This truncation of the expansion is driven by the respective values of «2 and A*' ar SPS and RHIC, and is by 
no means mandatory. 
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Figure 1. Conditional probability distribution /(pt2|Pti) of the relative angle of an associated 
particle for a trigger oriented along (full line) or perpendicular to (dashed line) the reaction plane. 



to find an associated particle close or back-to-back to the trigger (A9912 ~ or vr) is smallest. 
These rough features, entirely dictated by elliptic flow, are illustrated in figure [H where the 
values pTi = 6 GeV/c, pT2 = 4 GeV/c, V2{pt2) = 0-2, V2 = 0.1, = (500 MeV/c)^ and 
A'^ = 8000 have been used. In further detail, for both in-plane or out-of-plane trigger particles, 
the probability for associated particles is higher away (A(/?i2 ~ vr) than close (A(^i2 ^ 0) to 
the trigger: this is the trademark of the back-to-back correlation ([T]) induced by momentum 
conservation. 

Now, the aim of a correlation study would be to investigate this finer structure. For that 
purpose, one would like to "remove" the anisotropic-flow-induced pattern. Given the general 
expression of the conditional probability distribution (resp. joint probability distribution) as a 
function of the two-particle cumulant [1], the most natural — although experimentally highly 
challenging — recipe would be to divide /(PT2IPT1) [resp. /(pti;PT2)] by the single-particle 
distribution /(PT2) [resp. by the product /(pti)/(pt2)]) so as to isolate the correlation!!] Yet 
the usual procedure is rather to subtract a flow-modulated background, whose normalization is 
for instance fixed by requiring that the yield vanish at its minimum ("ZYAM" [9]). Following 
this approach, the "fiow-subtracted distribution of associated particles" is pictured in figure [2] 
for both in- and out-of-plane trigger particles. As anticipated, this distribution is larger away 
from the trigger than close to it. Furthermore, the away-side probability is smaller when the 
trigger points perpendicular to the reaction plane than when it points along this was also 
not unexpected given the shape of the two-particle probability distribution However, the 
detailed structures are quite non-trivial, in particular the dip at A(^i2 = vr in the case of an 
out-of-plane trigger, which both is a remnant of the "incompletely" subtracted anisotropic fiow 
patternll 

^ This procedure does not entirely suppress the influence of the anisotropic expansion, since the latter afi^ects 
the strength of the correlation according to the pair azimuth, see V2 in equation (O. However, inasmuch as it is 
experimentally feasible, it would indeed cancel the single-particle modulation induced by anisotropic flow, while 
leaving intact the azimuthal dependence of the two-particle cumulant. 

^ If /(PT2) had been divided from the conditional probability, instead of being subtracted from it, the resulting 
quotient would have been a smooth first-harmonic sinusoid, rather than the curve in figure (2] 
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Figure 2. "Flow-subtracted distribution of associated particles" /(PT2IPT1) ~ct/(PT2) '^s. the 
relative angle, for a trigger particle along (full line) or perpendicular to (dashed line) the reaction 
plane, a is fixed by the "ZYAM" condition. 

The azimuthally-dependent behaviours described above follow solely from assuming 1. that 
the particle transverse momenta sum up to and 2. that the particle emission is anisotropic: 
(Px) > (Py)- Those mere two ingredients are enough to give rise to non-trivial structures, 
which then have to be disentangled from those arising from additional physical effects. Two 
such patterns, whose analogueqil have reportedly been observed at RHIC, were identified: (a) 
a lower away-side {Aipi2 ~ vr) probability of associated particles in the case of an out-of-plane 
trigger compared to an in-plane trigger [lOj, and (b) a dip in the away-side yield of associated 
particles |1H [T2] . These or similar structures have been predicted in models of parton energy loss 
that involve some path-length dependence [13] [pattern (a)], or as signature of the interaction 
("Mach cone", "gluon Bresstrahlung" , "Cerenkov ring", "jet deflection", see e.g. reference [12j ) 
between the away-side parton and the medium through which it propagate [structure (b)]. If 
such models are to yield quantitative results, the possible contribution to the data of "trivial" 
correlations induced by momentum-conservation has to be investigated seriously. 
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